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Abstract 

Wc discuss the phase transitions between three states of a plasma fluid (plasma ball, uniform 
plasma tube, and non- uniform plasma tube), which are dual to the corresponding finite energy black 
objects (black hole, uniform black string, and non-uniform black string) localized in an asymptoti- 
cally locally AdS space. Adopting the equation of state for the fluid obtained by the Scherk-Schwarz 
compactification of a conformal field theory, we obtain axisymmctric static equilibrium states of 
the plasma fiuid and draw the phase diagrams with their thermodynamical quantities. By use of 
the fluid/gravity correspondence, we predict the phase diagrams of the AdS black holes and strings 
on the gravity side. The thermodynamic phase diagrams of the AdS black holes and strings show 
many similarities to those of the black hole-black string system in a Kaluza-Klein vacuum. For 
instance, the critical dimension for the smooth transition from the uniform to non-uniform strings 
is the same as that in the Kaluza-Klein vacuum in the canonical ensemble. The analysis in this 
paper may provide a holographic understanding of the relation between the Rayleigh-Plateau and 
Gregory-Lafiamme instabilities via the fiuid/gravity correspondence. 
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1 Introduction 

The AdS/CFT correspondence provides one of the best ways to explore the nonperturbative regimes 
of a certain class of strongly coupled quantum field theories by working in classical gravity. Conversely, 
from the viewpoint of black hole physics, we may be able to understand strong gravitational phenom- 
ena by studying an appropriate limit of field theories. In this paper, we investigate a class of AdS 
black holes, in particular, their phase structures and possible phase transitions, via the fluid/gravity 
correspondence (see, e.g., [1,2]), which has been proposed as a generalization of the AdS/CFT corre- 
spondence. 

Recently, Aharony, Minwalla, and Wiseman [3] argued that in a class of large-A gauge theories, a 
plasma ball (i.e., a lump of 'gluon' plasma) appears via a first-order deconfinement phase transition 
above a critical temperature. They also claimed that such a plasma ball will map to a finite energy 
black hole localized in the IR (infrared) region, and discussed what happens on the gravity side by 
studying various phenomena on the fluid side such as a plasma ball production and its subsequent 
decay by hadronization. In particular, they explicitly constructed a numerical domain wall solution 
to the Einstein equations with a negative cosmological constant which smoothly interpolates two 
spacetimes corresponding to the confined and deconfined phases. As a generalization of the analysis 
in [3], Lahiri and Minwalla [4] constructed a rotating plasma ring and a ball, which are dual to a 
rotating black ring and a rotating black hole, respectively (see also [5]). By construction (which we 
will briefly review in the text), these gravity duals reside in the IR region of AdS space compactified 
on a Scherk-Schwarz circle. This kind of black ring and hole solutions have not been discovered yet. 
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However, by comparing the phase diagrams of a plasma ring and plasma ball with those of a black 
ring and black hole in an asymptotically flat space, they found qualitative agreements between them. 

In this paper, we investigate the phase structures of another important class of black objects, 
i.e., black strings and black holes in a spacetime in which some spatial dimensions are compactified, 
from the viewpoint of the fluid/gravity correspondence (see [6] for the analysis on a similar system 
from a holographic viewpoint). We solve the d-dimensional relativistic Navier-Stokes equations [2] 
with an appropriate surface term to obtain static axisymmetric fluid equilibrium configurations. The 
configurations contain a spherical ball, a uniform tube, and a non-uniform tube. To see what their 
gravity duals are, we recall the arguments in [3]. We consider the following solution to the {d + 2)- 
dimensional Einstein equations with a negative cosmological constant, Rab = —{d + l)£~'^gab- 



= ^ e^"[-dt^ + T2n{u)de' + dwf] + T^^du' , (1) 
where i = 1,2, . . . ,d — 1, and £ [0, 27r). The function Tx{u) is defined by 

T,{u) = l-[-{d+l)e-\ . (2) 

This spacetime is regarded as an AdSd+2 space with a Scherk-Schwarz compactification, which is called 
the AdS soliton [7]. By taking u — > oo, one recovers the AdS(i4-2 in the Poincare coordinates with a 
uniform circle of 6. It should be noted that there is a cutoff in the IR region (i.e., small-n region) 
since the Scherk-Schwarz circle shrinks to a point at a finite value of u. Imposing a periodicity in the 
imaginary time, i.e., t = t + (3, {t = it), one can regard this spacetime as a thermal gas of gravitons 
at temperature T = (3~^ . In addition, there exists another exact solution 

dsl+2 = f (^e^^[-Tpiu)dt^ + de^ + dwf] + j^^du^) > (3) 

which has the same asymptotics as the spacetime (1). This solution can be regarded as a black brane 
with temperature T = I3~^. Considering both spacetimes, Eqs. (1) and (3), together in an ensemble, 
one can show that in the low temperature regime (T < Tc = l/27r) the AdS soliton (1) has a lower 
free energy, while in the high temperature regime (T > Tc) the black brane (3) has a lower free energy 
and dominates. Thus, the system undergoes a Hawking-Page type phase transition at T = Tc [7,8]. 

In [3,4], the plasma balls and plasma rings which extend uniformly in the ^-direction and non- 
uniformly in the u)j-directions (which are assumed to be noncompact) were constructed on the UV 
(ultraviolet) boundary by solving the fluid equations^. In this paper, we compactify one of Wj's on 
a circle, which we denote by coordinate z, and consider the plasma tube and plasma ball, where the 
former wraps the circle of the z-direction. As described in [4], the horizon topology of a dual black 
object is obtained by fibering the plasma configuration with an that shrinks to zero size at the 
fluid edges^. Therefore, the plasma ball B'^'~^ and plasma tube B'^~^ x S^, which we will obtain in 
this paper, map to the black objects whose horizon topologies are 5^ and S'^~^ x S^, respectively. 
Namely, they are a black hole and a black string, respectively, in an asymptotically AdS^ space on the 
Scherk-Schwarz circle Sg and Kaluza-Klein circle Sz- 



^Note that the dual black holes and black rings also extend in the M-direction, although they are localized in this 
direction. 

■^This condition on the shrinking of the circle corresponds to the fact that at the IR wall mapped from the outside 
of the plasma, the Scherk-Schwarz circle shrinks to zero size. 
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The gravity dual of the plasma tube obtained in this paper is rather different from the black 
string solutions in an asymptotically locally AdS space which have been studied extensively in the 
literature [9], as well as from the black strings in the asymptotically locally flat Kaluza-Klein space. We 
will see, however, that the thermodynamic properties of our AdS black hole-black string system display 
various qualitative similarities to those of the black hole-black string system in the asymptotically 
locally flat Kaluza-Klein space. The phase transitions and their dimensional dependence will be 
elaborated in this paper. 

This paper is organized as follows. In Sec. 2, we reduce the Navier-Stokes equations for a plasma 
fluid to a hydrostatic equation for axisymmetric equilibrium states. We then demonstrate that this 
equation is equivalent to that of Plateau's problem (or capillary minimizing problem), obtaining 
the constant mean curvature surfaces which model soap bubbles with a variational principle. In 
addition, we show that the uniform plasma tube is unstable against a perturbation, which is the fluid 
counterpart [10-12] of the Gregory-Laflamme instability [13]. In Sec. 3, an equation of state for the 
fluid is introduced which is obtained from the Scherk-Schwarz compactification of a conformal field 
theory. The thermodynamic variables for each phase are then calculated. In Sec. 4, the thermodynamic 
phase structures are obtained. The final section is devoted to a summary and discussion. A technical 
issue about numerical integrations is treated in Appendix A. A justification of numerical results and 
a natural interpretation of phase diagram are described in Appendix B. 

Note Added 

During the preparation of this paper, we were informed that a similar work had been completed 
independently [14]. It is interesting to compare their results with ours. 

2 Hydrostatics of Axisymmetric Plasma Lumps 
2.1 Relativistic Navier-Stokes Equations 

The equation of motion for the boundary fiuid is simply given by the conservation of the stress 
tensor [2], 

The stress tensor consists of the perfect fluid part, the dissipative part, and the surface contribution 
in a certain long wavelength limit as 

perfect surface dissipative' ' 

where the perfect fluid part takes the usual form as 

^pTrfect = iP + P>'^' + Pr (6) 

with the pressure P, proper energy density /o, and velocity fleld u^. The surface contribution can be 
written as 



TTurface = ^(^^^'^ " gn^/d¥^ (7) 

with a being a tension of the boundary. We have assumed that the surface of the fluid is given by 
$(x'^) = and the unit normal of this surface is denoted by := d^^{d^ ■ 9<1>)~^/^. For static 



3 




Figure 1: An axisymmetric static equilibrium state of fluid in a d 
schematically embedded in a three dimensional space. 



(n + 3)-dimensional flat spacetime 



fluids which we will study in this paper one can show the dissipative term does not contribute to the 
equation of motion. 

Now, we consider a, d = {n + 3)-dimensional flat spacetime (n > 1)^ in cylindrical coordinates 
= {t,z,r,(j)'), {i = l,2,--- ,n), 



-di^ + dz^ + dr^ + r'^-fij 
71abdx''dx' + r\,{(l))d(l)'d^, 



(8) 



where x'^ = {t,z,r), 7]ab = diag.(— 1, 1, 1) and 'jijdcp'^dcj)^ is the line element of the unit n-sphere. 
Non-zero components of the Christoffel symbol are 



-pa 



jk 



(9) 



where the bar is used for the variables with respect to the metric jij. Assuming that the fluid is static 
in this frame {u^ = 6^), the perfect fluid part of the stress tensor and its divergence are given by 



perfect 



/ p 
P 
OOP 

\ r-2pyj J 



perfect 



dP dP , 

oz or 



(10) 



where we have used an identity Vi-jjk = 0. The boundary surface of an axisymmetric fluid configura- 
tion can be given by ^{z, r) := r — h{z) = 0, where h{z) is a height function (see Fig. 1) assumed to 
be a single-valued function. Then, noting that the unit normal vector n'^ is given by 
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the surface contribution to the stress tensor and its divergence are given by 





(11) 



j^tiu ^ (^^{r - h) 

surface ^/Y^ 



/ 1 + 



V 



-1 
-h' 





-h' 





-r-2(l + /i'2)yi ) 



V^TTurface = ^'^(r - /i) ( 0, -(n + l)h' H , {n + l)H, 



(12) 



^In fact, the conformal boundary, where the plasma fluid resides, is [d + l)-dimensional if we take into account the 
existence of the Scherk-Schwarz circle, in which direction we assume the fluid configuration is uniform. Therefore, the 
total spacetime dimension of the AdS bulk is d + 2. 
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Here, H is the mean curvature of the surface defined by 



I ( h" n \ 

The first and second terms in the parenthesis correspond to the axial and azimuthal principal curva- 
tures of the surface, respectively (see, e.g., [12]). 
The non-trivial components of Eq. (4) are 

BP 

— - (n + l)a5{r - h)h' H = 0, 
dP 

— + (n+l)a5(r-h)H = 0. (14) 
or 

Away from the boundary r = h{z), we have the equation in the bulk, which gives 

P{z,r) = const. (15) 

Integrating the equations of motion across the boundary, we have 

P> - P< + (n + l)aH{z) = 0, (16) 

where P> and P< are the pressures just outside and inside the boundary, respectively. Now, we 
consider the case in which the pressure outside vanishes (P> = 0) for simplicity. Then, Eq. (16) reads 



n 



+ 1 V (1 + /i'2)3/2 hVT+h^J in + l)a 



This is the governing equation that determines the axisymmetric equilibrium states of the fluid. This 
equation implies that the fluid surface is given by a constant mean curvature surface. As we will see, 
this equation is derived by a variational principle which minimizes the surface area of the fluid while 
keeping the volume fixed. A uniform tube found as a trivial solution of Eq. (17) will suffer from the 
Rayleigh-Plateau instability^ . 

2.2 The Equivalent Plateau Problem and Rayleigh-Plateau Instability 

In general, the interior volume and surface area of a fluid body are written as 



V= dV, A= dA. (18) 

J interior J surface 

When we assume the staticity and axisymmetry of the fluid of which surface is given by r = h{z), the 
induced metric on the surface is given by 

dSsurfaee = "dt' + (1 + h'^)dz^ + /l^^,^- (^)d</.Mc/.^' . (19) 



*The axisymmetric constant mean curvature surfaces in general dimensions were obtained in [11], and their geometric 
properties were investigated to compare those of black strings. 
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Now, we assume that the z-direction is compactified on a circle, z E [— L/2,L/2]^. Noting that the 
volume and surface elements in Eq. (18) are given by dV^ = r"-y/7dr A dz A d^^ A • • • A dc/)"" and 
d^ = /i"v^(l + /i'2)7dz A d0^ A • • • A d(/)", we have 

»L/2 />L/2 



V\h\=^n\ dz/i'^+i(z), = (n + / dzVl+7^/i"(z), 

J-L/2 J-L/2 

^-^=— t/ V7d(A^A---Ad,/."= ^.w^^-ir (20) 
n + iy^n"^ r[(n + l)/2 + lJ 



The equation of motion for the configuration that minimizes (or extremalizes) the surface area for 
a given volume is obtained by varying the action 

m:=a^A\h\-v^V\h\. (21) 

Here, the constant uo corresponds to the surface tension, while the constant po corresponds to the 
pressure, but it is just a Lagrange multiplier to mathematically ensure that the volume is fixed. 
The Euler-Lagrange equation obtained by variation (bhl = 0) is equivalent to Eq. (17) with the 
identification of po ^ P< and ao ^ a. 

It is noted that by adding a term corresponding to the rotation energy to the action (21), which 
is roughly given by —ljJ with an angular velocity lv and an angular momentum J, the equation of 
h{z) for a rigidly rotating non-relativistic fluid is obtained. In the rotational case, the surface deviates 
from a constant mean curvature surface [4,5,10] ^ 

Equation (17) has two trivial solutions representing a uniform tube (UT) and a spherical ball (SB), 

h = hjjT := ro, H = i^uT : 



(n + l)ro ■ 



h = hsB := V^o-^^^ H = Hsb:=^. (22) 

As is well known, a translationally invariant cylindrical body is unstable if the linear dimension is 
longer than its circumference in d = 4 (the Rayleigh-Plateau instability). Its onset mode in general 
dimensions is obtained by the following static perturbation of the uniform tube [10, 11]. First, we 
expand h{z) around the uniform tube, 

h{z) =ro + ehi{z) + 0{e^). (23) 



""This compactification corresponds to compactifying a spatial direction of the field theory coordinates in the Poincare 
patch from the AdS point of view. However, it does not introduce the problem of singularities in the asymptotic IR region 
since the IR is cutoff by the horizon or the shrinking of the Scherk-Schwarz circle. We would like to thank T. Wiseman 
for discussion on this point. 

®See Ref. [14] for more general relations between the equation of motion (16), called the Young-Laplace equation, 
and variational principles. They showed that in a fully covariant way the Young-Laplace equation can be obtained by 
variational principle of either the entropy maximization or potential energy minimization. It was also shown that these 
two variational principles reduce to the area minimization with volume fixing in the static situation. Its special case 
corresponds to ours. Their general proofs nicely relate the Young-Laplace equation (derived from the Navier-Stokes 
equation), thermodynamic equilibrium condition (the entropy maximization), mechanical equilibrium condition (the 
potential energy minimization), geometrical condition (the area minimization). More importantly it also supports the 
original argument in [3] (and our thermodynamic consideration in Sec. 4), that the thermodynamics of finite energy 
black holes in the Scherk-Schwarz compactified AdS is mapped holographically to that of deconfined plasma lumps held 
by surface tension. The authors thank an anonymous referee for letting us know this point. 
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Substituting this expansion into Eq. (17), we have the linear perturbation equation at 0(e), 



h'l + ^ ^ ^hi = 0. (24) 

n 

With a boundary condition, say h'^ (0) = 0, we have 

hi{z) = /i(^) cos (knpz) , fcRP := (25) 

where h^^^ is an integration constant. Solution (25) corresponds to the marginally stable mode of the 
Rayleigh-Plateau instability. The uniform tube is unstable if the length of the cylinder L satisfies 
L > Lrp := 27r//cRp. In other words, the uniform tube is unstable if the radius rg satisfies tq < 
^'RP := y/nL /2'7r for a given L. This dimensional dependence of the critical mode is quite similar to 
that of the Gregory-Laflamme instability. See Refs. [10-12] for more on the similarities between the 
Rayleigh-Plateau and Gregory-Laflamme instabilities. 



3 Equation of State and Thermodynamic Variables 

In order to investigate the thermodynamic properties of the fluid lumps in the context of AdS/CFT, 
we introduce an equation of state for the plasma fluid. Then, the thermodynamic variables for each 
phase of fluid lumps are calculated. 

3.1 Equation of State 

The equation of state for the plasma is obtained by the Scherk-Schwarz compactiflcation of a (d + 1)- 
dimensional conformal fleld theory [3]. The free energy in terms of temperature T and volume V is 
given by 

J = (po - aT'^+i) V, (26) 

where po is a vacuum energy and a is some constant. The pressure P, entropy S, and energy £ are 
given by 

£ := 3" + TS = (^po + (iaT'^+^) V. (27) 
We then introduce the energy density p and entropy density s, given as 

P-= ^ = po + daT'^+\ 

, := I = (d + l)ai/('^+i) (po + P)'^/('^+^) . (28) 
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From the parameters po, a, and a, a characteristic length scale (whose existence violates the 
conformal invariance), temperature, and entropy density can be defined as 



^0 := — , 
Po 



poy/id+i) 

a J 



so 



(29) 



The temperature in Eq. (29) is nothing but the critical temperature of the confine-deconfine phase 
transition which occurs at / := 3"/V = —P = 0. 

Now, we consider the thermodynamic variables for the plasma lumps, that is, not only those for 
their constituent plasma but also the contributions from the surface term. Their energy density and 
entropy density are given by 



s = {d + l)so 



Tr 



(30) 



The energy, entropy, and Helmholtz free energy for the plasma lumps can be obtained by integrating 
the above densities, 



E = I (T^erfcct + T^surfacc) dV = pV + aA, 

S = j sdV = sV, 

F = E-TS = -PV + a A, 



(31) 



where we have used that the energy density and entropy density are constant in the present system. 

Prom the relations of P = —po + aT"""*"^ and P = {n + l)aH (note that d = n + 3), we obtain a 
relation between the temperature T and mean curvature H, 



T = T, 



1 + (n + l)loHL 



l/{n+4) 



J ._ ^0 



(32) 



We will discuss later a physically reasonable choice of the dimensionless parameter Iq. 

It is convenient to define the following normalized dimensionless variables when we draw phase 
diagrams, 

E ^ S F ^ T 



E--=-E^, S:=^, F:=—, T := — , (33) 
-C/RP JSB -TSB -tRP 

where E^ip and Trp are the energy and the temperature of the critical uniform tube, and Ssb and 
-FsB are the entropy and free energy of the spherical ball, respectively, which are given in the following 
subsection. 



3.2 Uniform Tube and Spherical Ball Phases 

For the uniform tube, /iut(-2) = tq with the length L, we can obtain the thermodynamic variables as 
functions of single parameter ro/L, 



EvT = Po 



0„L 



n+2 



FvT = Po- 



n 



{n' +An + 1)T5+4 - 1 



n 



-'UT 



T. 



T\jT = Tc 



1 + nln 



l/(n+4) 



(34) 
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where Tut := Tut/Tc- The critical values for the RP instability are obtained just by setting tq = trp 
(= y^L/27r) in Eq. (34). For instance, the energy and temperature of the critical tube are given by 



Erp = po 



n 



2tt 



n+l 



(n2 + 4n + l)flp 



1 



, Trp = Tc(i + lo 



l/(n+4) 



(35) 



where Trp := Trp/Tc. 

For the spherical ball, /isb(-z) = \/ ~ -^"^ ™ period L, we have the thermodynamic variables 
as functions of a single parameter Rq/L, 

n+2 



EsB = Po 



n + l 



(n^ + 5n + 5)T<: 



■n+4 
SB 



1 



-5*33 = so{n + 4)0„+iT 



^SB I ^ 



n+2 



^SB = /OO 



n+l 



^SB 



Rq 
L 



l + (n+l)iol|;) 



-1 l/(n+4) 



(36) 



where TgB := Tsb/T^. 

According to the definition (33), we have the following dimensionless quantities, 

27^y+^ (n2 + 4n + l)[l + n[o(T/ro)]-l fro\n+i 



S\JT 

Tut ■ 



a 



n+l 



(n2 +4n+ 1)(1 + 27rV^Zo: 

r ,r , ^ 1 ("+3)/("+4) 
1 + 7ilo{L/ro 



l + (n+l)Zo(T/i?o 



1 VT 



n"0. 



r?in+4 



1 



(n + l)"+il]„+i 
for the uniform tube phase, and 
nOn+i 



^0 



Tut 



l + n[o(T/ro) 
1 + 2TTy/n Iq 



n+2 

RoJ 

l/(n+4) 



(37) 



EsB 



l + (n+l)[o(T/i?o 



n2 + 5n + 5)[1 + (n + l)lo{L/Ro)] - I ( Rq 

L 



1 + 27r-v/n /o 



(n2 + 4n + 1)(1 + 27r^ Zq) - 1 

l/(n+4) 

, 'S'SB = TsB = 1 



n+2 



(38) 



for the spherical ball phase. 

Now, for each value of Iq, which should be specified, we obtain thermodynamic relations such as 
S = S{E) and F = T(T) implicitly via the parameter tq/L ^ [0, +oo) for the uniform tube phase and 
via Rq/L G [0, 1/2] for the spherical ball phase, respectively'^. Note that due to the fact Rq < L/2 
(i.e., the condition that a ball fits the period), the spherical ball phase has an upper bound of energy 
and a lower bound of the temperature. 



^Because we normalize the entropy by that of the spherical ball, Eq. (33) , the expression of Sut in Eq. (37) contains 
the radius of the spherical ball _Ro- When one draws the {E\jt,Svt) curve by varying the parameter ro, this Rq should 
be regarded as a function of r^. Such a relation _Ro ~ Roiro) is obtained by the combination of _Ro ~ Ro{E) obtained 
by inversing the relation in Eq. (38) and E — E{ro) in Eq. (37). Similar transformations are needed to draw the 
(-BnuTiSnut) curve [see S'nut in Eq. (43)]. 
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3.3 Non-uniform Tube Phase 

As given in the preceding sections, Eq. (17) is solved by the uniform tube and the spherical ball. In 
addition to these two solutions, Eq. (17) allows the third solution, which is a non-trivial constant 
mean curvature surface corresponding to the non- uniform tube (NUT) [11]. For the d = 4 case, this 
solution is known as Delaunay's unduloid [15], which is the surface of revolution of an elliptic catenary 
[see Fig. 2(a)]. Although a certain parametric representation of the curve is known for the d = 4 
case, for d > 4, we have to integrate Eq. (17) numerically. As we will see below, if we introduce 
a 'non-uniformness parameter' A := r„/r+, where r__ and are the smallest and largest radii of 
a non-uniform tube, all non-uniform tube solutions are described by a one-parameter family of A, 
(0 < A < 1). 

First, we introduce the following dimensionless variables, 

y := Hqz, w := Hoh{z), (39) 

where Hq is the mean curvature of the non-uniform tube. Then, the first integral of Eq. (17) is written 
in a potential form, 

(^pUi.)=0, UH:=l-[^^)\ (40) 

where K is an integration constant. We assume that there are two zero points of the potential U{w) 
denoted by w;+ and W- (0 < < w^). w+ and w- correspond to the maximum and minimum 
(dimensionless) radii of the non-uniform tube, respectively. Thus, A = W-/wj^ and A ^ 1 
correspond to the ball and the critical tube, respectively [see Fig. 2(b) for a schematic picture of 
U^w)]. If we first give a A, from U{w±) = 0, we obtain 




Figure 2: (a) The non-uniform tube in = 4 is called Delaunay's unduloid in geometry, which is the 
surface of revolution of an elliptic catenary. The elliptic catenary is obtained by rotating an ellipse along a 
line (i.e., the z-axis) and tracing the focus, marked by small solid circles, (b) A schematic picture of the 
potential U{w), defined in Eq. (40), for three values of A = w^/w^. From the top curve to the bottom, 
the value of A decreases within the range < A < 1. 
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per the period in terms of A and Hq. 

2 




zrn+l 



V{X). 



(42) 



As a result, for the non-uniform tube, we can obtain the thermodynamic dimensionless quantities 
as functions of the non-uniformness parameter A, 



-E^NUT 

5nut ■ 
Fnut ■ 
Tnut '- 



n 



2tt 



n+l 



[jn + 3)V + i] [1 + (n + l)!oL] + V-A 
[(n2 + 4n + 1)(1 + 2tt^/E Iq) - l]L"+2 ' 

{n+3)/(n+4) ^ ^ ^ 



1 + (n + l)/oL 
5^n+i [l + {n + l)lo{L/Ro) 

n + 1 (rj^+^ - inA - V) 



V 



- ,-Y 



^^ri ( rpn+4 



n+2 ^ _ y 



(43) 



where Tnut := ^nut/^c, and L{X), A{X), and F(A) are the dimensionless variables defined in Eq. (42). 
For each value of Iq, which should be specified, we have thermodynamic relations, 5 = S{E) and 
F = F{T), implicitly via the non-uniformness parameter A S (0, 1), (see also Footnote 7). 

Now, we have expressed all thermodynamic variables as functions of the parameter A in the in- 
tegration form. Although the integrations in Eq. (42) are of course finite, the integrands diverge at 
the both ends of the integration range since U{w±) = 0. Therefore, one needs some manipulations 
for an accurate numerical integration, especially in the highly deformed regime A <C 1. We present 
some technical prescriptions making it possible to figure the fine structures of phase diagrams in 
Appendix A. 



4 Phase Structure 



In the previous section, we have described the necessary thermodynamic variables of each phase by 
the appropriate parameters: the uniform tube is parameterized by tq/L G [0, oo); the spherical ball by 
Rq/L G [0, 1/2]; and the non-uniform tube by A G (0, 1). To draw phase diagrams, we have to specify 
the undetermined parameter Iq = Iq/ L = a/{poL). Although we are treating the sharp boundary as a 
fluid surface, in reality the boundary surface has a thickness of order a / po = Iq [3,4]. Therefore, 

we obviously should work in the limit Iq <^ L, where the thickness of the boundary can be ignored 
(except for highly deformed configurations). Otherwise higher-derivative contributions to the surface 
stress tensor would dominate. Hence, we take Iq = 1.0 x 10"^ in all concrete examples hereafter. It 
should be stressed, however, that the qualitative aspects of the phase diagrams do not depend on a 
specific choices of parameter Iq. In particular, we have confirmed that the critical dimensions which 
we will find do not change for other choice of /q within the range of 1.0 X 10-2 < [q < 1- It is noted 
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that such a choice of small Iq corresponds to < Trp/Tc — 1 <C 1 as one can see from Eq. (35). That 
is, we only consider the plasma just above the critical temperature Tc. 

4.1 Microcanonical Ensemble: S = S{E) 

First we discuss the phase transitions among the plasma lumps in the microcanonical ensemble. Hence, 
we consider the E-S relation. According to the behaviors of each phase curve, e.g., the number of 
cusps appearing in the E-S diagram here, we divide the dimensions into three groups, i.e., 4 < d < 9, 
10 < d < 12, and d > 13. We now discuss the possible phase transitions between the three phases 
(a spherical ball, and uniform and non-uniform tubes) with the criterion that the maximum entropy 
state for a fixed energy is favored. 

4.1.1 4<d<9 

In this first class of dimensions, the branches of the uniform tube and spherical ball intersect at the 
point labeled by A in Fig. 3. The spherical ball always has a larger entropy than the critical tube. 

From Fig. 3, one can see that the branch of non-uniform tubes emerges from the Rayleigh-Plateau 
critical point, located at -E = 1, and reaches the end of the spherical ball branch located at the point 
E, where the ball 'touches itself at the boundary z = ibL/2. The most important point is that the 
non-uniform tube branch always has a smaller entropy than the other two. Thus, the phase diagram 
suggests that the non-uniform branch is never favored. Suppose that we have a fat uniform tube (i.e., 
large E) and decrease its energy E. As E decreases, the uniform branch meets the spherical ball 
branch at the point A, and at this point the uniform tube transits to a ball with a discontinuous jump 
in configuration. This transition is of first order. 

Here, we can see that the energy-entropy diagram in d = 5 is quite similar to the mass-entropy 
digram of the black hole-black string system in the asymptotically locally flat 5 and 6 dimensional 
Kaluza-Klein space (see Fig. 6 in [16]). Only one apparent difference is that the localized phase, i.e., 
the spherical plasma ball does not deform unlike a localized caged black hole. 



d = 5 
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A/ SB „ 
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0.999 


VUT 








12 3 4 

A 


5 


6 



E 



Figure 3: The energy-entropy diagram for d = 5, containing the phases of the spherical ball (SB, blue- 
solid line), uniform tube (UT, black-dashed curve), and non-uniform tube (NUT, red-solid curve with small 
circles of data). The maximum entropy state for a given energy is favored. There is no cusp on the 
non-uniform tube phase, which is the case for 4 < d < 9. 
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4.1.2 10<d<12 



In this class of dimensions, the phase diagrams are somewhat comphcated due to the appearance of 
two cusps in the {E, S) diagram. The {E, S) diagrams for d = 10, 11 and 12 are shown in Fig. 4. The 
uniform tube and spherical ball intersect in d = 10 at the point A, while they do not intersect any 
more in d = 11 and 12 (in fact, ^d > 11). The spherical ball has a larger entropy than the critical 
tube in d = 10 and 11, while this is not the case in d = 12 (in fact, > 12). As can be seen from 
these observations, there appear rich and non-trivial phase structures in these transient dimensions. 

From the fact that the non-uniform branch has two cusps in the {E, S) diagrams, the non-uniform 
tube branch can be divided into three branches (i.e., three smooth curves): the first branch leaves the 
Rayleigh-Plateau critical point; the second one is connected to the end point E of the spherical ball 
branch; the third one is the curve between the first and the second ones. The first and third branches 
correspond to the non-uniform tube branch appearing in the lower dimensions (4 < d < 9). We shall 
first discuss the case of d = 11. The non- uniform tube emerges from the Rayleigh-Plateau critical 
point. E increases at first, however, it turns to decrease, and then to increase again. Eventually the 
non-uniform tube branch reaches the end point E of the spherical ball branch. During this change, 
the non-uniform tube branch intersects with the uniform tube branch at the point C and with the 
spherical ball branch at the point B. Thus, if we begin with a fat uniform tube (i.e., large E) and 
decreases E, the uniform tube transits to the non- uniform tube at the point C with a discrete jump in 
configuration. As E decreases further, the non-uniform tube transits to a spherical ball at the point 
B with a discrete jump in configuration. 

The phase structures of d = 12 is quite similar to that in d = 11. That is, if we begin with a 
fat uniform tube and decrease the energy, the favored phases shift as (uniform tube) — > (non-uniform 
tube) — > (spherical ball). These two phase transitions are accompanied by the discrete jumps in 
configuration, and hence they are of first order. However, the phase structure in d = 10 is different 
from those of d = 11 and d = 12, although there are two cusps in all cases. The branches of uniform 
tube and spherical ball still intersects (at the point A) . The spherical ball has a larger entropy than 
the critical uniform tube. As a consequence, the behavior of the non-uniform tube branch is similar 
to those in 4 < d < 9 except that two quite small cusps appear at the middle of the branch. That is, if 
we begin with a fat uniform tube and decreases the energy, the uniform tube transits to the spherical 
ball at the point A. The non-uniform tube has nothing to do with this transition. 



(a)d=10 (b)d=ll (c)d=12 




1 1.5 0.8 0.9 1 0.5 1 

A A A 

E EE 

Figure 4: The energy-entropy diagrams for (a) d = 10, (b) d = 11, and (c) d = 12. There appear two 
cusps on the non-uniform tube phase in this class of dimensions. 
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d=13 




0.5 1 



Figure 5: The energy-entropy diagram for d = 13. There exists one cusp on the non-uniform tube phase, 
which is the case for d > 13. 



4.1.3 d>13 

As mentioned before, the branches of the uniform tube and spherical ball do not intersect each other in 
this class of dimensions. The spherical ball phase has a smaller entropy than the critical uniform tube 
in this class of dimensions. The cusp near the Rayleigh-Plateau critical point disappears at d = 13, 
while the cusp near the spherical ball phase remains (we have confirmed this up to c? = 15). Therefore, 
if we begin with a fat uniform tube (i.e., large E) and decrease E, the non- uniform branch meets the 
Rayleigh-Plateau critical point, and at this point it smoothly transits to the non-uniform tube, which 
has a larger entropy. This transition is not accompanied by a discrete jump in configuration, and 
therefore is of second or higher order. As E decreases further, the non-uniform branch intersects with 
the ball phase at the point B, and transits to a ball. This transition is of first order. 

Thus, the smooth transition from the uniform tube to the non-uniform tube is realized for d > 13. 
We conclude that the critical dimension in the microcanonical ensemble is cf^i^rocan _ -^^ Thus, the 
critical dimension for the gravity dual is expected to be around /)™icrocan _ ^^microcan _|_ 2 = 

4.2 Canonical Ensemble: F = F{f) 

In the canonical ensemble, according to the number of cusps appearing in the (T, F) diagram, we again 
divide the dimensions into three groups {4 < d < 9, d = 10, and d > 11). The temperature-free energy 
diagrams for d = 5, 10, and 11 are shown in Fig. 6. Note that as mentioned before, the spherical 
ball phase has a lower bound of temperature, which is realized when the ball touches itself at the 
boundary. One can see how the transitions between the three phases occur by repeating discussions 
similar to those in the previous subsection. 

In d = 5, if we begin with a fat uniform tube (corresponding to a small T) and increase T, the 
uniform tube transits to a spherical ball at the point A, which is a first-order transition. The non- 
uniform tube is not involved in this transition. This transition pattern in the canonical ensemble is 
similar to that in the Kaluza-Klein black hole-black string system of low dimensions^. 

*Here, we adopt the definition of critical dimension by Sorkin [17]. However, if we define a critical dimension, from a 
non-perturbative viewpoint, by the dimension at and above which the transition from a uniform tube (or uniform black 
string) to a non-uniform tube (or non-uniform black string) is realized, we find 1)'^'"°'^'^'^ — 13, 

®We could not find literature in which the temperature-free energy diagram of the black hole-black string system is 
presented based on actual numerical data. But, see Fig. 1 in [6] for a temperature-free energy diagram in 10 dimensions 
speculated from a perturbation analysis, which exhibits a similar behavior to Fig. 6(a). 
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Figure 6: The temperature-free energy diagrams for (a) d = 5, (b) d = 10, and (c) d = 11, containing 
the phases of spherical ball (SB, blue-solid line), uniform tube (UT, black-dashed curve), and non-uniform 
tube (NUT, red-solid curve with small circles of data). The phase having the smallest free energy is favored 
for a given temperature. There is no cusp for d = 5 (in fact, 4 < < 9), while there exist two cusps for 
d = 10. One cusp is found in d = 11 (in fact, d > 11). 



The (T, F) curve of the non-uniform tube in d = 10 has two cusps. However, the non-uniform 
tube branch is always above the uniform tube branch. Therefore, the non-uniform tube branch has 
nothing to do with a reahstic transition, as the d = 10 case in the microcanonical ensemble. That is, 
only the first-order transition from a uniform tube to a spherical ball occurs (at the point A). 

In d = 11, if we begin with a fat uniform tube and increase the temperature, the uniform tube 
meets the Rayleigh-Plateau critical point, located at T = 1, and at this point it transits smoothly to 
the non-uniform tube, which has a smaller free energy than the uniform tube phase. This transition 
is of second or higher order. As the temperature increases further, the non-uniform tube meets a 
spherical ball at the point B. At this point, it transits to a spherical ball via a first-order transition. 

In the canonical ensemble, the critical dimension at and above which the smooth transition from 
the uniform to non- uniform tube phases is realized is d'^^'^ = 11. Hence, the critical dimension of the 
dual black strings is expected to be around D^^^ = d'^'^^ -|- 2 = 13. 

The patterns of the phase transitions in both microcanonical and canonical ensembles are summa- 
rized in Table 1. 



Table 1: Diverse patterns of the phase transitions among the uniform tube (UT), non-uniform tube (NUT), 
and spherical ball (SB) in the microcanonical and canonical ensembles. The '2nd' means that the order of 
transition may be either second or higher. 



Spacetime Dimension on Fluid Side: d — n + 3 


4 - 9 


10 


11 


12 


13- ... 




Transition Type 


UT — 


SB 


UT — > NUT — > SB 


Microcanonical 


Order of Transition 


1st 


1st & 


i 1st 


'2nd' & 1st 




No. of Cusps in {E, S) 





2 


1 




Transition Type 


UT — 


SB 


UT — > NUT — > SB 


Canonical 


Order of Transition 


1st 


'2nd' & 1st 




No. of Cusps in (T, F) 





2 


1 


Spacetime Dimension on Gravity Side: D — d + 2 


6 - 11 


12 


13 


14 


15 - ... 
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5 Summary and Discussion 



Adopting the equation of state for the fluid obtained by the Scherk-Schwarz compactification of {d + 
1) -dimensional conformal field theory, we have investigated the thermodynamic properties of the 
'deconfined gluon plasma' lumps, which are expected to be dual to black holes and strings localized 
in the IR of AdS^ on the Scherk-Schwarz and Kaluza-Klein circles. We have invoked the fluid/gravity 
correspondence in order to predict the phase diagrams of the AdS black holes and strings. We have 
found that those phase diagrams are qualitatively similar to those of the black hole-black string system 
in the asymptotically locally flat Kaluza-Klein space. These results are not surprising in the sense that 
it had been known that the phase structure of the black hole-black ring system in AdSs (AdSg) with 
Scherk-Schwarz compactification is qualitatively similar to that of the black hole-black ring system 
in the 5(6)-dimensional flat background [4,5]. However, there is as yet no clear explanation for the 
agreement between the phase structures of such distinct systems. Furthermore, the critical dimensions 
found in this paper ^/jmicrocan _ and -DJ^"^ = 13) are 'very close' (indeed equal in the canonical 
ensemble) to the ones in the asymptotically locally flat Kaluza-Klein space {D'^^^'^^ = 14 and 
-^♦KK ~ ^'^) respective ensembles [17,18]. This closeness/coincidence, respectively, may stem 

from a universality of critical dimension, or there may exist other unknown reasons. As mentioned 
in the Introduction, the bulk duals may have a non-trivial dependence on the holographic radial 
coordinate u. Furthermore, the distribution of the size of the Scherk-Schwarz circle on horizon cannot 
be determined a priori before solving the Einstein equations. However, the n-dependence as well as 
the distribution of ^-circles may be simply uniform in the limit of the large black holes/strings for 
some reason, and in this case the phase structures would resemble those of the Kaluza-Klein system. 
It would be interesting to clarify the reason for the qualitative agreement, by solving the Einstein 
equations to find the gravity duals explicitly. 

Next, we comment on the validity of our calculation, in particular, under what conditions the 
effective fluid description of the field theory is valid. One condition is that the length scale over which 
the temperature and pressure vary must be larger enough than the mean free path of quasiparticles 
of the system, which is of the same order as the inverse of the deconfinement temperature, (in 
the large 't Hooft coupling limit) [1,4]. In our system, the temperature and pressure are constant 
throughout the configurations, and this condition holds. Another condition is that the temperature 
must not be far from the critical temperature. Otherwise, it is no longer valid to assume that the 
surface tension is a constant, whose value was implicitly assumed to be that at the critical temperature 
a = (t{Tc) in this paper. From Fig. 6, we find that the temperature of non-uniform tubes always 
satisfies Tnut/Trp = 0{1). Since we have set Trp > Tc (by taking *C 1), Tnut/^c = 0(1) holds 
for all configurations near the transition points we are interested in. Finally, the radius of curvature 
of the fiuid surface in all directions must be much larger than the thickness of the surface, which is of 
order T~^. If this does not hold, higher-derivative contributions to the surface stress tensor must be 
included. Since we have discussed the case L ^ Iq = a/po ^ T~^, we can safely ignore the thickness 
of the surface. However, the waist of non-uniform tube can become arbitrarily thin in the limit 
A = r_/r+ 0, where the thickness of the surface cannot be ignored. Thus, the transition process 
from a non-uniform tube to a spherical ball could be modified by higher-derivative contributions. One 
may wonder whether this problem can be resolved by a field theoretic approach. 

One may object that the critical dimensions found in this paper are too large to occur in a 
realistic setting such as M/String theory, and cannot describe reality. However, a critical dimension 
at and above which a stable non-uniform phase appears will depend on various factors of the system 
under consideration. For instance, the critical dimension for boosted black strings is expected to 
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be significantly lower for highly-boosted strings [19]. Taking into account the fact that a large black 
ring [20] is regarded as a boosted black string, it is expected that there is a stable undulating black ring 
phase even in low dimensions such as 5 and 6. Therefore, it would be quite interesting to investigate 
the existence of the Rayleigh-Plateau instability for rotating plasma rings [4, 5] as well as its possible 
nonlinear consequences^'^. In a similar vein, it would be also interesting to investigate the Rayleigh- 
Plateau instability of ultra-spinning plasma balls (see [22] for a related observation on the gravity 
side), and to examine the (in)stability of the non-uniform tube phases obtained in this paper. These 
analyses may provide holographic interpretations of the stability of higher dimensional black holes. 
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A Manipulations for Accurate Numerics 

The integrands in Eq. (42) diverge since U{w±) = 0. Therefore, one needs some manipulations to 
numerically integrate them accurately, which we describe here. 

The variables in Eq. (42), L{X), A{X), and V{X), can be written as 

X(A)= Pdu;4==> {X = L,A,ovV), (44) 

Jw- yJ-U{w) 



where ipxiw) is a certain regular function, i.e., Tpiiw) := 2, tpAiw) := \J\ — U {w), or ipvii^) ■= 
2w^~^^. Using the explicit form of the potential U{w), one can write the integral as 

Jw- y (u)" -|- + K){w-^ — w){w — W-)g{w) 
where we have defined (n — l)-st order polynomial g{w) by 

„ ,n+l _ n _|_ Tf 

gH = 9.^w- := ^^-^y (46) 

m=0 ^ ' ' ^ ' 



^"See [21] for an interesting discussion on the (non-) existence of Saturn like configuration of plasma and its implications 
on the gravity side. We would like to thank one of the authors, C. Krishnan, for informing us of his paper. 
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Comparing the coefficients of botli sides of Eq. (46), we have a recursion relation of 

Qm - {w+ + w^)gm+i + = 0, (0 < m < n - 3), 

K {w+ + w-)K ^ ^ 

50 = , 91 = —7 gn-2 =W++W--1, gn-1 = 1- (47) 

This recursion relation is solved as 

g^ = 3n-2 - _ 9n-2 - ^+gn-l ^n-l^n^^ (2 < m < n - 3). (48) 

tt;+ — W- — W- 

Using the explicit expressions of K and w± , and Eq. (47) , we obtain the final expression, 

Y^ylTT) ' (0<m<n-l). (49) 

Then, we change the integration variable by 

w = + (50) 

Thus, the integral (44) is written as 

xw = t de (£i±£lM^^, (51) 

7-1 ^(u;" + 'u;'^+i +K){1- e)g{w) 

where w is regarded as the function of by Eq. (50), and g{w) is given as the polynomial through 
Eqs. (46) and (49). Finally, we use the following expressions for {K,'w±) to reduce the round-off error, 
which are equivalent to the original definitions (41), 



K;+ = ^r°. , ^_ = ^=°^^ , / . (52) 

e:;=oA"^ e™=oA'" (e::.=oA-)"^' 

We find that the expressions (51) and (52) significantly reduce the numerical errors and make it 
possible to obtain non-uniform tube solutions, even for A ^ 1. A mathematica notebook is available 
at the author's website, http://www.phys.huji.ac.il/~umpei/. 



B Justification and Interpretation of Cusp Structures 

Here, we argue that the cusp structures observed in the phase diagrams are not numerical artifacts. The 
phase curves were obtained by numerical integration of Eq. (51), which was done with NIntegrate 
in the mathematica code. We confirmed the convergence of integrated values by changing the 



Table 2: Numerical values of the slenderness parameter, Eq. (53). 



d = n + 3 
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6 
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10 


11 


12 


13 


14 


A 


3.142 


2.221 


1.814 


1.571 


1.405 


1.283 


1.187 


1.111 


1.047 


0.9935 


0.9472 
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parameters of AccuracyGoal, PrecisionGoal, and WorkingPrecision. That is, the digits required to 
draw the fine structures of phase diagrams such as the cusp structure do not change against the change 
of these parameters. It is noted that such a convergence is made rapid by virtue of the manipulations 
described in Appendix A. The fact that the curve of non-uniform tube phase approaches the Rayleigh- 
Plateau critical point and the tip of the spherical ball phase at the ends of curve supports our numerics 
independently. 

We also note that the existence of cusp structures, especially during the transient dimensions 
discussed in Sees. 4.1.2 and 4.2, is not inconsistent with our conclusion that the phase structure of 
black strings and black holes in the Scherk-Schwarz compactified AdS is similar to that of the black 
hole-black string system in the asymptotically locally flat Kaluza-Klein space. That is, the critical 
dimensions known in the Kaluza-Klein black hole-black string system [17, 18] were predicted by higher- 
order perturbations, which tells us the behaviors of a phase curve only near the Gregory-Laflamme 
critical point. Therefore, although the critical dimension of Kaluza-Klein black string has been known 
to be -D*,KK = 14 in the microcanonical for example, the possibility is not excluded that the non- 
uniform black string branches in -Drk = 12 and Z^kk = 13 behave like the curves of non-uniform tube 
in Fig. 4 (note that the fully non-linear behaviors of non-uniform black string phase curves have been 
obtained for Dkk < 11 in [23], but such a behavior was not observed). Therefore, the cusp structures 
in the transient dimensions found in this paper are not only consistent with the known gravitational 
calculations but also suggest an important lesson that a critical dimension should be determined non- 
periurbatively, or lesson that the critical dimension defined by the smooth transition from the uniform 
branch to non-uniform branch is not enough to know the global structure of a phase diagram. It would 
be interesting to perform numerical analyses to obtain the fully non-linear behaviors of non-uniform 
black strings for -Dkk = 12 and 13, although their difficulty was one of the original motivations to 
work on fluids in this paper. 

In Ref. [11], we found the critical dimension at and above which the non-uniform tube branch 
emanates from the Rayleigh-Plateau critical point with decreasing volume V for a fixed period L, of 
which existence is related to the critical dimension in the thermodynamic phase diagrams obtained 
in this paper. We found a simple criterion whether the non-uniform branch increases or decrease its 
volume near the Rayleigh-Plateau critical point as follows. Let us define a 'slenderness' parameter of 
the critical uniform tube. 



We may say that a critical tube is slender, say, if A > 1, while one is fat if A < 1. The numerical 
values of this parameter are given in Table. 2. From the table, the critical tube is found to be slender 
for 4 < d < 12, while fat for d > 13. On the other hand, the non-uniform tube branch increases its 
volume for 4 < < 12, while decreases it for d > 13 near the critical point (see Figs. 2 and 4 in [11]). 
The above agreement of two threshold dimensions seems not to be just a coincidence from a simple 
geometric consideration: the non-uniform tube branch emanating from a slender critical tube has to 
increase its volume for a fixed period L in order to approach a spherical ball, while one emanating from 
a fat critical tube has to decrease its volume. Thus, the slenderness provides us a simple criterion 
to forecast the behavior of non-uniform tube around the Rayleigh-Plateau critical point, and our 
numerical results support this geometric expectation and fit our intuition. 

Finally, it should be stressed that the transient behaviors of non-uniform tube observed in Fig. 4 
(and also Fig. 6) are quite natural in that the drastic change of phase structure between the lower 
dimensional class (4 < d < 9) and higher dimensional class {d > 13) was 'smoothly interpolated' by 



A : 



2ro 




TT 



(53) 



19 



the appearance of the 'swallowtail' (i.e., the part of the curve between the pair of two cusps) at d = 10 
and its subsequent growth at d = 11, 12. It would be interesting to regard dimension d as a fictitious 
continuous parameter and perform the same calculations to obtain phase diagrams for non-integer 
'dimensions' 9 < d < 13, which would show us a continuous change of phase structure. 
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